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Preface 

These  notes  contain  the  material  presented  during  a 
series  of  lectures  at  the  Massachusetts  Institute  of  Technology 
in  May,  1975.   It  was  the  intent  of  these  lectures  to  famil- 
iarize the  students  with  the  basic  concepts  of  MHD  stability 
theory  and  to  lead  them  to  a  point  where  they  are  able  to 
understand  and  interpret  the  literature  on  this  subject. 
Therefore,  only  a  rather  limited  number  of  applications  is 
given  v/hich  are  by  no  means  representative  for  the  present 
status  of  the  theory. 

Due  to  its  relative  simplicit:y ,  MHD  theory  can  be  used 
to  study  the  effect  of  txhe  geometry  of  the  magnetic  field 
upon  plasma  confinement.   Since  this  question  is  crucial  for 
the  problem  of  Controlled  Thermonuclear  Fusion,  MHD  is  one 
of  the  most  extensively  studied  subjects  in  plasma  theory, 
and  there  exists  a  vast  literature.   Nevertheless,  it  should 
be  emphasized  that  MHD  stability  is  not  very  well  understood 
at  present.   For  one  thing,  it  is  not  clear  how  one  should 
define  stability  in  order  to  obtain  a  problem  which  is  both 
physically  relevant  and  mathematically  tractable.   Moreover, 
it  has  not  been  possible  to  gain  a  generally  valid  qualita- 
tive picture  within  any  of  the  common  definitions  of 
stability.   Correspondingly,  it  is  not  exaggerating  to  say 
that  the  relevant  questions  cannot  be  answered  without 


extensive  numerical  computations.   The  practically  important 
knowledge  thus  consists  largely  of  more  or  less  unrelated 
special  cases  which,  nevertheless,  are  extremely  involved. 
Naturally,  none  of  these  can  be  presented  here. 

Due  to  the  somewhat  unusual  point  of  view  taken  in  these 
lectures,  much  of  the  material  included  has  been  gathered 
from  relatively  inaccessible  or  even  unpublished  sources. 
Therefore,  and  consistently  with  the  attem.pt  to  keep  these 
notes  self-contained,  a  list  of  references  is  not  included. 
It  is  hoped  that  the  numerous  original  authors  accept  this 
apology. 
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1.   Some  Simple  Consequences  of  the  Nonlinear  Equations 


We  write  the  nonlinear  equations  of  ideal  MHD  in  the  form 

du 

Prrf  =-Vp   +   J    X   B  (1.1) 

dt  '^        ^        ^ 


J   =   curl   B  (1.2) 


div   B   =   0  (1.3) 

^  =   curl    (u   X   B)  (l.'^) 

at  ~       ~ 

1^  +   div    pu   =    0  (1.5) 

o  t  — 


^     ^     =   0  (1.6) 

dt      pY 


Here  3/3t  is  the  time  derivative  at  a  fixed  position,  and 
d/dt  is  the  time  derivative  following  the  fluid  motion: 


^      -      ^  +  u-V  (1.7) 


dt     9 


1/ 
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These  equations  are  not  sufficiently  tractable  to  serve  as 
a  basis  for  a  stability  theory.   Therefore,  we  consider  their 
linearized  versions  in  most  of  what  follows.   However,  we  derive 
a  few  general  consequences  which  help  us  to  develop  some  intu- 
ition. 

Conservation  of  Energy 

The  equations  (1.1-7)  imply  the  identity 


^(^■P  u^  +  :^P  +  I  B^)  (1.8) 


1^  ..2  ,  _Y_ 

Y-1' 


+  divK^-P  u   +  v^p)u  +  B  X  (u>=B)]  =  0 


Integrating  this  over  a  fixed  spatial  domain,  we  caxi 
use  Gauss'  theorem  to  show  that  the  second  term  does  not  contri- 
bute if  we  adopt  the  boundary  conditions 


u   =  0   ,   B   =  0  (1.9) 

n  n 


(which  are  appropriate  for  a  rigid  and  perfectly  conducting  wall). 
A  subscript  n  denotes  the  normal  com.ponent  of  the  vector.  V/e  thus 
have 


d_ 
dt 


dT,(|pu^  +  ~jp  +  I"  B^)  =0  (1.10) 


The  three  terms  in  the  integrand  are  the  kinetic,  thermal, 
and  magnetic  energy  densities.   The  integral  of  the  first  is  the 
kinetic  energy,  (denoted  by  K)  ,  and  the  integral  of  the  sum  of  the 
of  the  last  tv,'0  is  the  potential  energy  (denoted  by  U)  . 


Frozen-ln  Magnetic  Field  Lines 

We  show  that  magnetic  field  lines  can  be  viewed  as  moving 
with  the  fluid  velocity  u.   We  use  the  fact  that  equation  (1.3) 
is  equivalent  to 

B  =  Va  X  V3.  (1.11) 

It  is  trivially  seen  that  equation  (1.11)  implies  equation  (1.3) 
but  to  show  the  converse,  viz.,  that  any  field  satisfying 
equation  (1.3)  can  be  represented  by  two  "flux  functions" 
a  and  3  according  to  equation  (1.11),  requires  some  thought. 
Equation  (1.11)  implies  that  both  a  and  3  are  constant  along 
field  lines.   Hence  each  field  line  can  be  viev\fed  as  labeled 
by  a  pair  of  values  of  a  and  3.   If  a  and  3  depend  on  time, 
each  field  line  moves  accordingly.   To  prove  the  assertion 
we  show  that  equation  (1.11)  implies  equation  (1.^)  if  we 
let  a   and  3  depend  on  time  according  to 


^  =  0    ,    If  =  0.  (1.12) 

dt  dt 


We  have 


^   =   |^(VaxV3)  =  Vff  X  V3+Vaxv|l 


=  -7(u-Va)xV3-VaxV(u-V3)  . 


On  the  other  hand 

curl(uxB)  =  curl[ (ux(VaxVe)] 

=  curl[(u-V3)Va-(u.Va)V6] 

=  V(u.V3)xVa-V(u-Va)xVB. 

Combining  of  the  last  two  relations  yields  equation  (1.4). 

The  statement  that  field  lines  move  with  the  fluid 
velocity  is  not  an  absolute  one  because  there  is  a  whole 
class  of  velocity  fields  u,  characterized  by 

curl[(u  -  u)xB  =  0, 

which  can  be  viewed  as  field  line  velocities.   The  point 
is  that  u  is  a  possible  field  line  velocity.   The  sifnificance 
of  this  is  that  an  identity  can  be  assigned  to  each  field 
line  for  all  times.   Therefore,  field  lines  can  bever  break 
or  combine.   This  is  often  referred  to  as  conservation  of  the 
topology  of  field  lines.   For  instance,  a  closed  field  line 
will  be  closed  forever,  two  closed  field  lines  will  be  linked 
in  the  same  way  forever,  and  field  lines  which  lie  in  some 
sui-face  will  lie  in  a  similar  surface  forever. 


circulation  Theorem 

The  circulation  along  some  closed  path  is  defined  as 

C  =  Idl'U  (1.12) 


In  isentropic  ideal  hydrodynamics  C  is  conserved  for  any 
closed  path  if  this  path  is  carried  with  the  fluid  motion. 
In  isentropic  ideal  MHD,  C  is  conserved  for  paths  which  are 
formed  by  magnetic  field  lines.   We  call  the  plasma 
isentropic  if  p  =  np   with  constant  n •   The  surfaces  of 
constant  pressure  then  coincide  with  those  of  constant 
density. 

To  prove  the  circulation  theorem  of  MHD,  let  us  first 

derive  an  expression  for  the  time  derivative  of  an  integral 

/ 
of  the  form  <)d£-A,  where  h-   is  any  time  dependent  vector 

field,  and  the  integral  is  along  a  closed  magnetic  field 

line. 


It  is  convenient  to  use  the  flux  functions  a  and  3, 
supplemented  by  a  third  function  x  subject  to  dx^dt  =  0, 
as  coordinates.   Then 


dZ   =  dx  B/J, 


(1.13) 


where 


J  =  (VaxVB)Vx  =  B-VX, 

Since  these  coordinates  are  moving  with  the  fluid,  the 
operator  d/dt  equals  the  time  derivative  at  fixed  a,^,    and  Xj 
and  thus  commutes  with  the  x-integraticn .   Hence 


A-B 


dx 


dt 


A-B 

J^dt^~  t^       J    dt^ 


Using  equations  (1.3-4)  and  the  relation 


dJ    ,  , . 
■7-;-  -   -J   div  u 
dt 


(1.14) 


(which  can  be  verified  by  a  straightforward  calculation) 
we  find 

dt^^§^-  -J-   dt 


dA 

=   B«  (~  -  uxcurl  A  +  V(A'U)) 
~   at    ~       -,      -,  ^ 


Now,  since  the  last  term  Integrates  to  zero. 


dt  ) 


r  SA 

d^-A  =  "d£  •  (tt^  -  ux  curl  A)  (1.15) 


For  the  circulation  theorem,  we  put  A  =  u  to  obtain 

8u 


dC    f     '^^ 


Using  Eq .  (1.1),  the  integrand  can  be  written  as 

^y  -  uxcurl  u  =  -V(-^  c^^h   ^  u^  )  +  i  J  X  B 

The  first  term  then  does  not  contribute  to  the  integral 

because  it  is  a  gradient,  and  the  second  term  does  not 

contribute  because  it  Is  perpendicular  to  the  path  of  in- 
tegration. 
Well-posedness 

Before  we  go  on  to  the  linear  theory,  we  would  like 
to  stress  one  important  point:   the  nonlinear  MKD  equations 

are  synmietric  hyperbolic,  and  therefore  the  initial  value 

problem  is  well-posed.   More  precisely,  for  any  initial 

data  which  are  compatible  with  the  equations,  and  with  the 

boundary  conditions  (1.9),  the  equations  possess  exactly 

one  solution   which  depends  continuously  on  the  initial 

data.   Even  though  we  do  not  work  with  the  nonlinear 

equations,  we  keep  this  fact  in  mind  whenever  we  encounter 

difficulties  with  the  equilibrium  equations  or  with  the 

linearized  equations:   one  might  sometimes  be  tempted  to 


conclude  that  the  MHD  model  Is  not  good  If  such  difficulties 
occur;  the  point  which  we  are  making  here  is  that  the  MHD 
model  is, at  least  from  the  mathematical  point  of  view,  as 
good  as  a  model  can  be,  and  difficulties  arise  only  if  we 
ask  unreasonable  questions.   It  is  reasonable  to  require 
that  the  temporal  evolution  of  the  system  is  uniquely 
(and  continuously)  determined  by  its  initial  state;  but 
there  is  no  logical  necessity  for  the  system  to  relax  to 
an  equilibrium  state,  nor  for  equilibrium  states  to  exist 
at  all. 


2  .  Static  Equllibi-'lum 

An  equilibrium  Is  a  time  Independent  solution  of  the 
equations  of  motion.  For  simplicity,  we  restrict  ourselves 
to  static  equillbriaj  i.e.  equilibria  with  u  =  0.   The 
equations  for  these  are 

Vp  =  JxB  (2.1) 

J  =  curl  B  (2.2) 

div  B  =  0  (2.3) 

An  immediate  consequence  of  equation  (2.1)  is 

B.Vp  =  0  (2.4) 

This  means  that  the  magnetic  field  lines  lie  in  pressure 
surfaces,  which  thus  are  "magnetic  surfaces".   For 
confined  equilibria,  these  surfaces  must  be  closed.   It 
can  be  shovm  that  equation  (2-3)  is  not  compatible  with 
closed  magnetic  surfaces  unless  each  surface  encloses  a 
multiply  connected  domain,  such  as  a  torus.   This  is  why  one 
is  primarily  interested  in  toroidal  equilibria.   In  a 
toroidal  equilibrium  the  pressure  surfaces  form  nested 
sets  of  toroids.   In  a  confined  equilibrium  the  outermost 
surface  carries  the  lowest  pressure.   It  is  convenient  to 
assume  that  this  outermost  surface  coincides  with  a 


perfectly  conducting  rigid  well.   The  magnetic  field  lines 
(and  also  the  field  lines  of  the  current  density)  lie  In 
the  toroidal  surfaces. 

It  can  be  shown  that  a  "rotation  number"  can  be  as- 
signed to  each  magnetic  field  line.   This  is  the  limit  of 
the  ratio  of  the  number  of  v;indings  the  short  and  the  long 
way  around  of  a  magnetic  field  line  as  one  follows  it 
forever.   The  rotation  number  is  the  same  for  all  field 
lines  in  the  same  surface.   The  field  lines  in  a  particular 
surface  are  closed  if  the  rotation  number  is  rational  there, 
and  each  field  line  covers  the  surface  densely  otherwise. 
For  instance,  the  field  lines  close  upon  themselves  after 
one  circuit  the  long  way  around  without  encircling  the  short 
way  if  the  rotation  number  is  zero;  but  they  close  the  short 
way  around  without  encircling  the  long  way  if  the  rotation 
number  is  infinite.   A  surface  with  rational  rotation  number 
is  called  "rational  surface."   If  the  rotation  number  varies 
from  one  surface  to  another,  the  equilibrium  is  said  to 
have  "shear".   The  rational  surfaces  are  then  dense,  but 
almost  all  surfaces  are  "irrational".   A  special  class  of 
shearless  equilibria  (the  rotation  number  is  constant)  are 
those  with  all  field  lines  closed  (the  rotation  number  is 
a  rational  constant).   Closed  line  equilibria  are  special 
in  many  respects,  and  their  theory  is  particularly  well 
developed. 

The  equilibrium  equations  (2.1-3)  are  neither  elliptic 
nor  hyperbolic:   the  magnetic  field  lines  are  (degenerate) 
characteristics,  but  there  is  also  an  imaginary  cone.   As 
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a  consequence,  the  theory  of  the  equilibrium  equations  (2.1-3) 
is  by  no  means  understood.   It  has  been  possible  to  formulate 
well  posed  boundary  value  problems  for  simply  connected  domains 
(the  real  characteristics  intersect  the  boundary),  for  toroidal 
domains  v/ith  axial  symmetry  (the  hyperbolic  part  of  the  equations 
can  be  integrated,  so  that  the  problem  reduces  to  a  purely 
elliptic  one),  and  for  certain  equilibria  with  zero  rotation 
number  (the  real  characteristics  are  closed  curves).   However, 
it  seems  likely  that  it  is  not  possible  to  formulate  a  well 
posed  boundary  value  problem  for  arbitrary  toroidal  domains 
because  of  the  complicated  topological  structure  of  the  real 
characteristics  (dense  set  of  rational  surfaces).   In  other 
words,  it  seems  likely  that  an  equilibrium  does  either  not 
exist  for  arbitrary  domains,  or,  if  it  happens  to  exist,  does 
not  depend  continuously  on  the  parameters  describing  the  shape 
of  the  domain.   In  the  present  notes  we  will  not  be  concerned 
with  equilibrium  theory,  but  rather  adopt  the  point  of  view 
that  some  equilibrium  is  given,  and  then  investigate  motions 
near  this  equilibrium. 
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3 .   Strict  Linear  Stability 

We  say  that  an  equilibrium  Is  "stable"  If  the  plasma, 
according  to  the  nonlinear  equations  of  motion,  stays  near 
this  equilibrium  forever  If  It  was  near  initially.   Unfortu- 
nately, this  definition  Is  neither  unique  nor  useful.   It 
Is  not  unique  because  there  are  many  ways  of  defining  what 
one  means  by  "near".   For  Instance,  one  could  require  that 
the  deviations  from  their  equilibrium  values  of  all  quantities 
appearing  in  the  equations  of  motion  are  small  everywhere, 
or  one  could  require  that  the  kinetic  energy  is  small.   The 
above  definition  is  not  useful  because  v;e  are  not  able  to 
handle  the  nonlinear  equations  well  enough  to  decide  whether 
a  given  equilibrium  is  stable  or  not. 

A  major  simplification  consists  in  linearizing  the  equations 
about 'the  equilibrium.   One  writes  B  =  3   +  B, ,  u  =  u  ,  etc., 
where  quantities  with  a  subscript  zero  satisfy  the  equilibrium 
equations,  and  then  neglects  terms  which  are  quadratic  in  the 
perturbing  quantities.   The  equations  (1.1-6)  then  become 

%  ^  -^Pl^i^^V-Io^^l  >  ^3.1) 

J^  =  curl  B^    ,  (3.2) 

div  B^  =  0   ,  (3.3) 

^=  curl(u^xB^)  ,  (3.i^) 

^Pl  +  div(PoUi)  =  0  ,  (3.5) 

at 

_.  =  -u^.Vp^  -  YP^  div  u^  .  (3.6) 
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This  is  a  set  of  linear  homogeneous  equations  for  the 
perturbing  quantities.   Like  their  nonlinear  counterparts, 
these  equations  define  a  well-posed  initial  value  problem. 
Nevertheless,  it  is  unknown  to  what  extent  the  solutions 
of  the  nonlinear  equations  are  approximated  by  solutions  of 
the  linearized  equations. 

Since  the  linearized  equations  are  homogeneous, 
smallness  of  a  solution  is  equivalent  to  boundedness.   We 
thus  call  an  equilibrium  stable  if  all  solutions  of  the 
linearized  equations  remain  bounded  as  t^-°°.   Of  course, 
we  have  to  decide  what  we  mean  by  "bounded".   In  other 
words,  we  have  to  introduce  some  norm  in  the  space  of 

vector  functions  (B,  ,u-,  , ),  and  then  call  a  solution 

bounded  if  this  norm  is  finite.   One  possibility  is  the 
maximum  norm.   We  then  call  the  equilibrium  stable  if  the 
quantities  B  ,u-,  , .  .  .  .  remain  everyv/here  bounded  as  t-*^ 
for  all  solutions  of  equations (3 • 1-6 ) .   To  distinguish 
this  definition  from  others,  we  call  it  "strict  linear 
stability."   This  definition  might  seem  reasonable  at 
first  glance,  but  unfortunately  virtually  no  equilibrium 
of  interest  is  stable  according  to  it. 

To  show  this,  we  pick  some  closed  field  line,  and 
define  the  first  order  circulation  as 

c(t)  =  f  d5,-u   .  (3.7) 

Since  the  equilibrium  is  static,  the  path  of  integration 
does  not  move  to  lowest  order,  and  the  time  derivative  of 
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c  is  Obtained  by  Just  differentiating  the  integrand.   Using 
the  linearized  equations  of  motion  and  the  equilibrium 
equations  (and  again  assuming  that  Vp^xVp^  =   q),  we  compute 


and  further 


dt  -  fl^'-Tt 


-  r^  fd£.(-Vp  +J  xB  +J  xB,  ) 
f^o  J  ~      1  -1  ~o  ~o  ~1^ 


dfc 
dt2 


If     1   3B, 

=   -  ^Jai   g^  Vp^.curl(u^xB^) 

=  -  ^JdJl  gl  div((u^xB^)xVp^) 

=  -  ^|d£gidiv(u^.Vp^)B^) 
=  -^{^.^/b^-V(u^.Vp^) 
=  -  ^  jd£.?(u^.Vp^)  =  0 


(3.8) 


(3.9) 
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Hence  dc/dt  Is  constant,  and 

c(t)  =  c(t^)+||(t-t^)  •  (3.10) 

The  initial  data  can  be  chosen  such  that  dc/dt  =j=  0  unless 
Vp   =  0  at  the  field  line  under  consideration.   Then  jcl->-«> 
and  hence  also  max  |  u,  | -»-«'  as  t-^=°j  and  the  equilibrium  is 
unstable.   The  only  equilibria  v;hich  cannot  be  shown  to  be 
unstable  in  this  manner  are  those  which  either  have  no  closed 
field  lines  at  all  or  in  which  the  pressure  gradient  is  zero 
at  each  closed  field  line  (which,  for  reasons  of  continuity, 
requires  that  the  pressure  is  constant  wherever  the  rotation 
number  is  not  an  irrational  constant). 

It  is  interesting  to  note  that  dc/dt  =  0  for  those 
perturbing  magnetic  fields  for  which  a  single-valued  vector 
field  4  can  be  found  such  that 

B^   =    curl(§xB^)  (3.11) 

For  these,  the  expression  (3-8)  for  dc/dt  is  identical  to  that 
for  d^c/dt^,  [eq.(3.9].  with  §  replacing  u.^.   If  eq.(3.11)  holds, 
then  C  can  be  interpreted  as  the  displacement  of  magnetic  field 
lines:   A  field  line  of  the  equilibrium  field  B^ ,  when  shifted 
by  an  amount  ^,    becomes  a  field  line  of  the  total  field 
B  =  B   +  B-,  .   To  demonstrate  this,  we  represent  the  equilibrium 
field  by 

~o      0    o    ' 
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and  then  shov;  that,  to  first  order  in  B, ,    the  total  field 
satisfies 

B  =  VaxV3, 

where 

a(x+|)  =  a^(x),3(x+p  =  3q(x), 

provided  equation  (3.11)  holds.   The  computation  proceeds  as 
follows : 

B  =  Va(x)xV3(x)  =  Va^(x-pxV3Q(x-|) 

=  V(a^(x)-C.Va^(x))xV(3o(x)-|-V3Q(x)) 

=  Va  xVB  -Va  xV(C-V3  )-V(5.Va  )xV3 

0     0     0-0      ~     0      0 

=  Va  x3^+  curl((C-V3)Va^-  (.^■Va^)V&^) 

=  Va  xV3  +  curK^xCVa^xVB^)) 

=  B   +  curl(?  X  B  ) . 
-0        -3    ~o 

The  fact  that  the  linear  circulation  c  is  a  linear 
function  of  time  seems  to  contradict  the  conservation  of  the 
nonlinear  circulation  C.   This  apparent  contradiction  is 
resolved  by  observing  that  c  does  not  approximate  C  in  general 
This  is  because  c  is  computed  by  integrating  along  a  field 
line  of  the  unperturbed  field  B  ,  while  C  is  computed  by 
Integrating  along  a  field  line  of  the  total  field  B   +  B,. 
These  two  paths  of  integration  are  unrelated  in  general  be- 
cause given  a  closed  field  line  of  B   there  need  not  be  a 

closed  field  line  of  B   +  B^  which  approximates  the  former. 

~o    ~1 

To  give  an  example,  we  consider  the  fields  B   =  V6,  and 
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and  B-,  =  ee,  where  G  is  the  angle  about  a  straight  axis,  and 

e  Is  the  unit  vector  along  this  axis.   Then  the  field  lines 

of  B  are  circles,  but  those  of  B   +   e  B,  are  helices  which 
~o  -o       ~1 

depart  arbitrarily  far  from  the  circles  even  for  arbitrarily 

small  G.   In  the  light  of  these  remarks,  the  fact  that  the 

linear  circulation  c  is  conserved  if  Eq .  (3-11)  holds  agrees 

with  the  fact  that  ^  then  is  a  displacement  vector  by 

means  of  which  a  relation  can  be  established  between  field 

lines  of  B  and  field  lines  of  B   and  B^  so  that  c  indeed 
~o  ~o     ~1 

approximates  C. 

What  can  we  now  conclude?   We  cannot  conclude  that  all 
plasmas  are  ill-behaved  because  experimental  evidence  shows 
the  contrary.   As  already  mentioned  in  Sec.  1,  we  can  also 
not  conclude  that  the  MHD  equations  are  no  good.   Thus, 
there  are  two  possibilities:  either  we  work  with  the  non- 
linear equations,  or  we  try  to  find  a  better  definition  of 
stability  within  the  linearized  MHD  model.   Obviously,  the 
second  possibility  is  the  only  feasible  one.   This  possibility 
is  further  Justified  by  two  facts:   Firstly,  confinement 
experiments  are  not  intended  to  last  forever;  therefore,  it 
is  crucial  to  find  out  how  much  plasma  and  hov/  much  energy 
are  lost  due  to  instabilities  v;ithin  the  intended  duration 
of  an  experiment.   Secondly,  the  MHD  equations  describe 
reality  poorly  in  certain  situations  (e.g.,  if  the  magnetic 
field  changes  appreciably  within  one  Larmor  radius);  therefore 
the  spatial  structure  of  an  Instability  has  to  be  known  in 
order  to  decide  whether  it  has  to  be  taken  seriously.   For 
these  reasons,  a  more  detailed  description  of  the  solutions 
of  the  linearized  MHD  equations  is  called  for. 
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^ .   stability  of  a  System  with  One  Degree  of  Freedom 

To  introduce  the  various  concepts  commonly  used  in  MHD 
stability  theory,  but  whose  significance  and  interrelations 
are  not  fully  understood,  we  consider  a  completely  tractable 
problem.   The  one-dimensional  motion  of  a  single  particle 
subject  to  a  conservative  force  Is  governed  by  the  equation 

X  =  F(x)  (4.1) 

Here  x(t)  is  the  position  of  the  particle  at  time  t,F(x)  is 
the  force  acting  upon  the  particle  at  the  position  x,  and  dots 
denote  time  derivatives.   The  motion  is  uniquely  determined 
for  all  times  if  both  the  position  x  and  the  velocity  x  are 
given  initially.   A  particular  position  x   is  an  equilibrium 

if  the  force  is  zero  there,  F(x^)  =  0.   The  equilibrium  is 
nonlinearly  stable    if  both  the  velocity  and  the  deviation 
from  the  equilibrium  position  remain  small  forever  provided 
they  are  small  initially,  or  more  precisely,  if 


|x(tQ)-x^i->0    and    \k{t^)\^0 


imply  that 


lx(t)-x  1^0  and  lx(t)j->0  for  all  t  >  t^ 


18 


In  terms  of  the  potential  energy  U(x)j  defined  by 

F(x)  =  -U'(x),  (4.2) 

an  equilibrium  is  characterized  by 


U'(x^)  =  0,  (4.3) 


and  stability  of  this  equilibrium  is  equivalent  to  the 

requirement  that  U  has  a  local  minimum  at  x  .   This  statement 
is  called  "nonlinear  energy  principle."   To  indicate  the  proof, 

we  note  that  equation  (4.1)  can  be  integrated  once  to  yield 

the  conservation  of  energy: 


|^(|x2+U)  =  0  (4.4) 


A  further  integration  then  yields  the  general  solution  in  the 
implicit  form 

X 

dx 
t-t   =      /2(E-U(x))  (4.5) 

1  2 
where  E  =  p"^  +U  is  the  (constant)  energy.   From  this  it  is 

seen  by  inspection  that  x(t)  and  x(t)  oscillate  about  x   and 

0  respectively,  with  amplitudes  going  to  zero  if  x(t  )  -»-  x 

and  x(t  )  -^  0 ,  if   U  has  a  minimum  at  x  ,  but  that  a  solution 

with  x(t  )  =  X   and  x(t  )  =  e  departs  away  from  x^  by  an 
o     o        o  o 

amount  0(1)  even  if  e->-0  otherwise. 

Next  we  discuss  the  linearized  theory.   When  linearized  about 
an  equilibrium  x  ,  the  equation  of  motion  (4.1)  becomes 

V  +  U"(x^)C  =  0  (4.6) 
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where 


^  =  ^-^o  (^.7) 


is  the  displacement  from  the  equilibrium.   As  before,  we  call 
the  equilibrium  "strictly  linearly  stable"  if  all  solutions  ^ 
are  bounded  as  t->«>.   A  more  optimistic  definition  is  that  of 
"exponential  stability":   no  solution  grows  exponentially  with 
time.   A  related,  and  still  more  optimistic  definition  is  that 
of  "a  stability":  no  solution  grows  faster  than  exp(ot).   To 
be  precise, 

lim  C(t)  exp(-ctt)  =  0 

for  a>o    in  a  CT  stable  equilibrium,  and  exponential  stability  is 
equivalent  to  a  stability  f or  a  =  o .   Obviously,  a^  stability 
Implies  o^    stability  if  o^    >   a^.      Therefore,  the  concept  of 
a  stability  is  more  flexible  than  that  of  strict  linear  stability 
(or  exponential  stability):   Rather  than  dividing  all  equilibria 
in  just  two  classes  labeled  by  "good"  and  "bad",  it  provides  a 
continuous  label  (viz.,  the  smallest  a  for  which  the  equilibrium 
is  a  stable) . 

With  the  ansatz  C°^  exp(itot),  equation  (^.6)  becomes 

(-aj^+U"(x^))5  =  0  (1].8) 

The  set  of  those  values  of  to   for  which  there  is  a  nontrivial 
solution  is  called  the  "spectrum"  of  the  system.   In  the  present 
simple  case,  the  spectrum  consists  of  the  single  real  eigenvalue 
w   =  U"(x^),  and  the  general  solution  of  equation  (4.6)  is 
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2 

If  to   =f  0  >  tiut 


C  =  a  exp(lwt)+b  expC-iwt)  (^.9) 


?  =  a+bt  (^.10) 


2 

in  the  exceptional  case  co   =0.   Here  a  and  b  are  complex 

constants.   To  be  physically  meaningful,  the  real  part  of  C 
has  to  be  taken. 

We  now  have  to  distinguish  three  cases: 

1)  U"  is  positive  at  x  •  Then  U  has  a  local  minimum,  the 
spectrum  is  positive,  and  all  solutions  of  the  linear- 
ized equation  of  motion  are  osciallatory .  The  equili- 
brium is  stable  according  to  any  of  the  above  definitions 

2)  U"  is  negative  at  x  •   Then  U  has  a  local  maximum,  the 
spectrum  contains  a  negative  eigenvalue,  and  there  are 
exponentially  growing  solutions.   The  equilibrium  is 
both  nonlinear ly  unstable  and  exponentially  unstable. 

It  is  a  stable  for  a  >_  /-U"(r  J,    but  a  unstable  otherwise. 

3)  U"  vanishes  at  x  .   Then  U  may  or  may  not  have  a  local 
minimum  depending  on  higher  derivatives,  the  spectrum 
is  non-negative  but  contains  the  origin,  and  there  are 
solutions  which  grow  linearly  though  not  exponentially. 

•The  equilibrium  is  strictly  linearly  unstable,  but 
exponentially  stable.   It  may  or  may  not  be  nonlinearly 
stable . 
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The  following  statements  sumrr.arize  what  we  have  learned 
so  far: 

1)  Exponential  stability  is  equivalent  to  U"(x  )  >0 
("linear  energy  principle")  and  a  stability  is 
equivalent  to  U"(x  )+o   >^  o  ("modified  energy 
principle") . 

2)  Strict  linear  stability  implies  nonlinear  stability, 
nonlinear  stability  implies  exponential  stability, 

and  exponential  stability  implies  o   stability. 

2 

3)  Strict  linear  stability  is  equivalent  to  oj  >0, 

2 
exponential  stability  is  equivalent  to  cj  >0,  and 

2   2 
o  stability  is  equivalent  to  to  +a  ^  0. 

To  give  a  preliminary  interpretation  of  the  fact  that 
most  MHD  equilibria  are  strictly  linearly  unstable  we  anti- 
cipate that  what  appears  to  be  the  exceptional  case  in  the 

2 
motion  of  a  single  particle  (viz.,cj   =  0)  is  the  common 

case  in  ideal  MHD.   We  also  anticipate  that  exponentially 

growing  MHD  instabilities  always  have  a  counterpart  in 

the  nonlinear  theory,  but  that  this  need  not  be  true  for 

non-exponential  instabilities. 
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5_^   Spectrum  and  Variational  Principles  In  MHD 

We  now  formulate  generalizations  to  ideal  MHD  of  the  concepts 
of  the  last  section.   This  is  tentative  to  some  extent  be- 
cause some  of  the  statements  which  are  trivially  true  for 
one  particle  have  to  be  postulated  for  a  plasma.   The 
reason  is,  of  course,  that  the  theory  of  partial  differen- 
tial equations  is  less  developed  than  that  of  ordinary 
differential  equations. 

Let  us  first  formulate  the  nonlinear  energy  principle 
of  ideal  MHD.   The  analogy  betv/een  the  energy  conservation 
laws  (1.10)  and  (4.^.)  suggests  to  postulate:   A  state 
is  an  equilibrium  if  the  potential  energy. 


U  = 


f  dt(-^p  +  ^  B^),  (5.1) 


is  stationery  along  all  possible  paths  of  the  system 
through  this  state;   an  equilibrium  is  "energetically 
stable"  if,  in  addition,  U  has  a  local  minimum  (again 
along  all  paths  through  the  equilibrium) .   A  restriction 
to  paths  of  the  system  was  unnecessary  in  our  mechanical 
problem  because  there  is  only  one  possible  path  in  a 
one-dimensional  motion.   If  we  omitted  this  restriction 
in  the  MHD  problem,  we  would  be  left  with  the  trivial 
equilibrium  p  =  0,  B  =  0.   It  is  exactly  this  restriction 
vjhich  makes  it  difficult,  if  not  impossible,  to  establish 
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the  relation  between  energetic  stability  and  some  (reasonably 

defined)  "nonlinear  stability".   The  reason  is  that  it 

does  not  seem  to  be  possible  to  describe  all  paths  through 

a  given  state  without  solving  the  equations  of  motion. 

A  convenient  way  around  this  difficulty  consists  in 

replacing  the  above  restriction  by  the  requirement  that 

certain  constants  of  the  motion  (viz.,  those  which  are 

known,  and  which  depend  only  on  p  and  B,  but  not  on  u) 

are  the  same  in  all  states  to  be  considered.   This  gives 

rise  to  a  different  definition  of  stability,  which  we 

call  "variational  stability":   An  equilibrium  is  variationally 

stable  if  it  minimizes  the  potential  among  all  states 

with  the  same  constants  of  motion.   This  statement  is 

called  "nonlinear  energy  principle". 

The  requirement  of  variational  stability  might  be 
more  restrictive  than  that  of  energetic  stability  because 
we  are  not  certain  to  know  all  relevant  constants  of  the 
motion,  so  that  we  might  admit  paths  v/hich  are  not  ac- 
cessible to  the  system.   On  the  other  hand,  in  the  in- 
stability exhibited  in  connection  with  the  circulation 
the  system  follows  a  path  which  does  not  pass  through 
the  equilibrium  (the  topology  of  the  field  lines  is 
different),  so  that  energetic  stability  might  be  an 
overly  optimistic  requirement.  In   summary,  variational 
stability  is  a  concept  whose  relation  to  the  solutions 
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of  the  equations  of  motion  Is  unknown.   This  is  so  despite 
the  fact  that  it  can  be  shown  that  the  variational  principle 
for  equilibrium,  using  constants  of  the  motion,  is  equivalent 
to  the  equilibrium  differential  equations  (2.1-3),  and 
despite  the  fact  that  there  is  a  relation  between  the 
linear  version  of  the  energy  principle  and  the  solutions 
of  the  linearized  equations  of  motion. 

The  nonlinear  energy  principle  is  tractable  to  some 
extent.   Correspondingly,  almost  all  the  work  v.'hich  ever 
has  been  done  in  nonlinear  MHD  stability  was  based  on 
the  nonlinear  energy  principle. 

Next  we  introduce  the  spectrum  of  linearized  MHD. 
As  before,  we  take  the  Fourier  transform  with  respect  to 
time  of  the  perturbing  quantities.   Time  derivatives  in 
the  linearized  system  (3.1-6)  are  then  replaced  by  a  factor 
i(jj.   As  a  consequence,  the  perturbing  density,  pressure, 
and  magnetic  field  can  be  expressed  in  terms  of  the  velocity 
u,  and  we  are  left  with  an  equation  for  u.   Omitting 
subscripts  we  can  write  this  equation  as 

w  pu  =  Pu  (5-2) 

where  F  is  the  matrix  differential  operator  given  by 

Fu  =  Bxcurlcurl(uxB)+  curl(uxB)xJ   -  V(u-Vp+YP  div  u),(5.3) 
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and  by  the  boundary  condition  u   =0.   The  perturbing  field 
B,  °^  curl(uxB)  then  Is  automatically  tangential  to  the 
boundary  If  the  equilibrium  field  has  this  property. 

It  should  be  noted  that  linear  growth  of  the  circulation 
now  has  been  eliminated  because  the  perturbing  magnetic  field 
automatically  satisfies  Eq.(3.11)  as  a  consequence  of  putting 
B,  °=  exp(lut).  The  motions  in  which  the  circulation  grows 
correspond  to  to  =  0,  the  only  case  in  which  elimination  of  B, 
and  P-,  is  a  restriction ,  and  which,  therefore,  is  not 
completely  covered  by  Eqs. (5.2-3). 

It  is  mathematically  convenient  to  imbed  the  solutions 

2 
of  equation  (5.2)  in  the  Hilbert  space  L  of  complex-valued 

vector  functions  u  which  satisfy  the  boundary  condition  and 

which  are  square  integrable, 

dT  I  u|   <  «>  . 

The  scaler  product  in  this  space  is  defined  as 

dx  u*«Up  , 


(u  ,u  )  =   , 
~±   -d  } }  J 


where  asterisks  denote  the  complex  conjugate.   The  spectrum 

2 
then  is  defined  as  the  set  of  those  values  of  ca   for  which 

2 
the  linear  operator  (F  -  u  p)  does  not  have  a  bounded 

Inverse.   According  to  the  alternative  theorem.,  those  values 

of  CO   for  which  equation  (5.2)  has  a  square-lntegrable 

solution  belong  to  the  spectrum.   These  are  called  point 

2 

eigenvalues.   However,  there  are,  in  general,  values  of  u 

which  belong  to  the  spectrum  but  are  not  point  eigenvalues. 
These  are  said  to  belong  to  the  continuous  spectrum.   If 
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2 

0)  belongs  to  the  continuous  spectrum,  a  sequence  u, 

can  be  found  for  which  llm(F-ca  p)u,  =  0,  but  which  itself 
does  not  converge. 

The  operator  F  has  the  property  of  being  "hermitian" , 
or  "formally  self  adjoint".   This  means  that 


(u^,Fu2)  =  (£;yi>U2^  ^^'^^ 

for  any  functions  u^  and  u„  in  the  domain  of  F.  (This 

2 

domain  is  smaller  than  L  because  not  all  square-integrable 

functions  are  twice  dif f erentiable ) •   The  property  (5-^) 
follows  from  the  identity 


~i'~~2~?~l'~2  (5-5) 


=  div[(u^xB)xcurl(u*xB)-(u*xB)xcurl(u^xB) 


)2    -'    ----^-1  f^  ^i^i 


'  c       ^ 


+YP(U2  div  u*-u«  div  U2)  +  (?' (f  i''!i2'*  ^^^ 

upon  integration,  and  using  Gauss'  theorem  to  show  that 
the  r.h.s.  vanishes  because  of  the  boundary  conditions. 
The  identity  (5.5)  can  be  verified  using  the  equilibrium 
equations  (2.1-3).   It  should  be  noted  that  the  operator 
F  is  not  selfadjoint  in  the  strict  sense  because  its 
adjoint  F* (defined  by  first  finding,  for  each  u,,  that 
u  for  which  (u-^^,  Fu^)  =  (u-.,u„)  identically  in  a^,    and 
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then  putting  u  =  F*u-,  )  has  a  wider  domain  than  F.   Never- 
theless, It  can  be  shown  that  F  has  a  self ad  joint 
extension.   This  is  crucial  because  the  pov/erful  theorems 
of  functional  analysis,  on  which  some  of  the  statements 
given  below  are  based,  would  not  apply  otherwise.   In 
what  follows  we  will  tacitly  assume  that  v/e  are  dealing 
with  that  extension  v/lthout  introducing  a  new  symbol  for 
this . 

As  a  consequence  of  the  self-adjointedness  of  F, 
the  spectrum  is  real.   As  we  will  see  later,  the  structure 
of  this  spectrum  is  quite  complicated:   There  are  both 
point  eigenvalues  and  continua,  and  the  former  frequently 
accumulate  at  the  tips  of  the  latter.   The  point  eigen- 
values are  governed  by  a  variational  principle.   In 
terms  of  the  two  quadratic  functlonals 

W[u]  =  |(u,  Fu)   ,   N[u]  =  |(u,pu)        (5.6) 

This  can  be  formulated  as 


w^  =  W/N   ,  6aj^  =  0,  (5.7) 


(where  the  symbol  6  denotes  the  variation  with  u )  j  or 
equlvalently 

W-w^N  =  0,   6(W-a)^N)  =0  (5-8) 
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The  first  version  of  this  variational  principle  can  be 
interpreted  as  requiring  that  W  be  stationery  if  u  is 
varied  subject  to  the  constraint  N  =  1.   The  second 
version  arises  from  observing  this  constraint  using 
a  Lagrangian  multiplier.      To  demonstrate  the  validity 
of  the  variational  principle  one  simply  shows  that  its 
Euler  equation  is  equatlon(5 • 2 ) • 

A  similar  variational  principle  holds  for  the  minimum 

2 

(denoted  by  w  .  )  of  the  spectrum: 
^         ^      mm 


0)^.   =  inf  S   >  C5.9) 

min       N 


or  equivalently 


co^.   =  max  I  oj^|W-aj^N>0  V  (5-10) 

min 


2 

If  u  .   is  a  point  eigenvalue,  this  is,  of  course,  a 
mm      ^       &       >        J 

special  case  of  the  principle  (5-7-8).   However,  it  happens 

2 

frequently  that  w  .   belongs  to  the  continuous  spectrum; 
^     ''       mm      ^ 

in  this  case  the  infimum  of  W/N  (i.e.,  its  largest  lower 

2     2 
bound)  is  not  attained  in  L  .   oa  .   can  be  shown  to  be 

'  mm 

finite  (i.e.,  >-°=)  as  long  as  both  |j|  and  Vp  are  bounded. 

On  the  other  hand,  u  .   is  never  positive  because  W[u]  =  0 

'   mm  ^ 

for  the  trial  functions  u  =  f(p)B+g(p)J  with  arbitrary  f 
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and  g.   Thus,  the  origin  is  always  a  point  eigenvalue  with 
infinite  multiplicity.   It  turns  out  that  it  also  is  the 

lower  edge  of  a  continuum. 

2 

We  now  call  an  equilibrium  "a  stable"  if  the  L   norm 

of  the  velocity  cannot  grov;  faster  than  expCcJt)  for  large 
t. 


lim||u||  exp(-at)  =  0  for  a>a  (5-11) 

t->-oo 


where  ||u||  =  / ( u , u ) .   Correspondingly,  we  call  it  "ex- 
ponentially stable"  if  it  is  a  stable  for  o  =  o .   This 
definition  is  natural  for  two  reasons:   Firstly,  it  relates 
to  the  spectrum,  and  secondly,  it  relates  to  the  kinetic 
energy  of  the  system.   As  to  the  relation  with  the  spectrum, 
it  can  be  shown  (using  Laplace  transform  techniques  and  the 

spectral  theorem  for  self-adjoint  operators)  that  o 

2       2 
stability  is  equivalent  to  u  ,   >_-  cr .   In  other  words, 

0   stability  obtains  if,  and  only  if,  the  spectrum  does  not 


extend  below  -a   (for  this  reason,  \/-w  .     is  called  the 

min 

"maximum  growth  rate").   The  variational  principle  (5.10) 
now  implies  that  ^  stability  is  equivalent  to 


W+0^N  >  0  (5.12) 
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This  is  the  "modified  energy  principle"  of  ideal  MHD. 

2 
For  a   =0,  this  reduces  to  the  ordinary  energy  principle, 

or  better  "linear  energy  principle":  Exponential  stability 

is  equivalent  to  W  >  0.   The  relation  to  the  kinetic 

energy,  defined  by 


K(t)  =  N[u(x,t)], 


arises  from  the  inequalities 


I  Pminl'y  il^    IK   <   -1-  p^^^   ||u  ||2 


These  imply  that  a  stability  is  equivalent  to 


11m  K  exp(-2at)  =  0  for  a   >    o  (5.13) 

t-^-oo 


We  have  not  proved  the  relation  between  the  spectrum  and 
o   stability,  as  we  have  not  proved  the  vaj^iational  principle 
for  the  maximum  growth  rate.   Instead,  we  will  give  an 
elementary  proof  of  the  equivalence  of  Eqs.  (5.12)  and 
(5.13)  which  is  independent  of  the  notion  of  the  spectrum. 

It  would  be  tempting  to  define  "strict  energetic 

stability"  as  positivity  of  the  spectrum,  i.e.,  as  co^ .   >  0. 

min 

However,  we  have  already  seen  that  w^ .   <_  0 .   Correspondingly, 
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it  seems  likely  that  perturbations  can  always  be  found 
for  which  K  -*■  «  as  t  -*■  <»  (in  some  non-exponential  way). 
This  has  not  been  shown  yet,  and  it  does  not  follow 
from  our  investigation  of  "strict  linear  stability" 
because  there  a  maximum  norm  was  used  as  opposed  to  the 
present  Hilbert  norm  (growth  in  the  Hilbert  norm  would 
imply  growth  in  the  maximum  norm,  but  the  converse  is 
not  true) . 

For  later  reference,  we  now  write  down  explicit 
expressions  for  the  functional  W.   A  simple  integration 
by  parts  yields 


W  =  I  ///dT||Q|2+(Jxu*)-Q  (5.14) 

+  (u'Vp)  div  u*+Yp|div  u|   (  , 


where 


Q  =  curl(uxB) 


A  particularly  useful  form  of  W  is 


1  C        J^'^p  ,? 

W  =  4  SffdT  \   |Q+ p  (u-Vp)r       (5.15) 


+  aIu-VpI^  +  YpI  div  u|^  J^   , 
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where 


2(JxVp) • (B-V)Vp 
A  = Z j^ .         (5.16) 

|Vp? 


For  real  u,  this  form  is  derived  from  Eq.(5-1^)  by  using 
the  identity 

(Jxu) • Q+(u-Vp)  div  u 

JxVp|2 
=  (u-Vp)^  (A+ 1^) 

+  2(u-Vp)(-?^^P)-g 
IVpl^ 


(u-Vp)(u-(JxVp)) 

+  div[  — = — ~ B]  , 

|Vp|2 


which,  in  turn,  can  be  verified  by  making  extensive  use 
of  the  equilibrium  equations. 
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6 .    Elementary  Proof  of  the  Modified  Energy  Principle 
We  wish  to  show  that 

11m  K  ep(-2at)  =  0  (6.1) 

t-voo 

for  a  >  Oj  and  for  all  solutions  of  the  linearized  MHD  equations 
(3.1-6),  is  equivalent  to 

V/^  =  W  +  a^N  >  0  (6.2) 

2 

for  all  test  functions  in  L  .   Our  proof  will  be  elementary 

in  that  no  reference  is  made  to  the  spectrum.   It  is  also 

noteworthy  that  it  will  suffice  to  assume  that  the  operator 

P  Is  hermltlan  (rather  than  strictly  selfadjolnt ) . 

We  relate  the  operator  F  to  the  linearized  equations 

(3-1-6)  in  two  distinct  ways,  without  taking  their  Fourier 

transform.   Introducing  a  "displacement  vector"  ^  by 

dK 
u  =  T^  (6.3) 

we  can  integrate  Eqs.  (3-4)  and  (3.6)  to  obtain 

B^  =  curl(CxB)+B  (6.^) 

P-L  =  -  C'Vp-YP  div  §  +  p,  (6.5) 

where  B(x)  and  p(x)  are  arbitrary  functions  (subject  to 
div  §  =  0),  and  subscripts  zero,  as  before,  have  been 
omitted  in  equilibrium  quantities.   Upon  substitution  of  the 
expressions  (6.4-5)  into  Eq .  (3.1),  this  becomes 


3^ 


3t2 


+  FC  =  -Vp  +(curl  B)xB+JxB        (6.6) 


Assuming  now  that  the  r.h.s.  of  Eq. (6.6)  is  zero,  we 
restrict  ourselves  to  the  sub-class  of  solutions  governed 
by  the  homogeneous  equation 

p   -  +  FC  =  0  .  (6.7) 

3t 

It  can  be  shown  that  these  solutions  are  exactly  those 
for  which  the  system,  passes  through  its  equilibrium  state. 

If  J  on  the  other  hand,  we  differentiate     Eq.(3.1) 
with  respect  to  time,  and  then  substitute  Eqs.(3.^)  and 
(3.6),  v?e  obtain 


B^u 


3t 


I  +  Fu  =  0  (6.8) 
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However,  only  those  solutions  of  Eq.(6.8)  are  associated 
with  solutions  of  the  original  system  (3-1-6)  for  which 
Eq.(3.1)  holds  initially ,  or  more  precisely,  for  which 
the  initial  acceleration 

a  =  p  3u/3t|^  ^  Q 
is  such  that 

-Vp  +  (curl  ExB+JxB  =  a 

can  be  solved  for  p  and  B. 
Loosely  speaking,  we  can  say  that  the  solutions  of 
Eq.  (6.8)  contain  those  of  the  actual  problem  (3-1-6), 
and  that  the  solutions  of  the  actual  problem  contain 
those  of  Eq.(6.7).   Therefore,  Eq.(6.7)  can  be  used  to 
demonstrate  necessity  of  the  modified  energy  principle 
(which  involves  exhibiting  one  solution  of  the  actual 
problem  whose  kinetic  energy  grows  like  exp(2a  )  if 
W   can  be  made  negative)  and  Eq.(6.8)  can  be  used  to 
demonstrate  sufficiency   (which  involves  showing  that 
the  kinetic  energy  grows  slower  than  exp(2at)  for  all 
solutions  of  the  actual  problem  if  W  ^0. 

Necessity 

Taking  the  scaler  product  of  Eq .  (6-7)  with  9?/3t 
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we  find 


|^(N[3|/at]  +  W[§])  =  0 


or  upon  integration, 

N[3§/9t]  +  W[C]  =  E.  (6.9) 

This  is  the  linearized  version  of  the  energy'  conservation 
(1.10).   If,  on  the  other  hand,  we  take  the  scaler  product 
of  Eq.(6.7)  with  C  ,    we  compute 


(C,p9^C/3t^)    =   |^(^,p9C/9t)-0§/3t,pa|/9t) 


T    d2 

i  S_p(CoP?)    -    OC/9t,p3?/9t) 

^  dt'^  ~   "" 

^2 

^  N[C]  -  2N[9C/9t] 
dt'^ 


to  obtain  the  "virial  equation" 

^-^   NEC]  -  2N[9C/9t]+2W[C]  =  0  (6.10) 

dt     ~  ~  - 

Upon  elimination  of  W  with  the  help  of  the  energy  conservation 
(6.9),  this  becomes 


^2 

2-,.  N[C]  -  MN[aC/9t]+2E  =  0  (6.11) 


dt 


We  now  assume  that  W  <  0  for  some   trial  function  ri,i.e., 
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I 


ee 


(51-9) 


^([5]N^)    <    [i]N^    [5    ]N 


u-fBCiqo    8M    (t7T"9)    P"^    (Sl'9)    ssTC^TiBnbauT    om:;    aqa   Suxuxqiuoo 


(t?T-9) 


I 


Xt^I/'  Md/)   ^(|d/id/)   >  ^(li  d/^5d/)   = 


se        5e 


5e 


z  ^p 


*s;n2   ^P- 


^P^ 


((5   d/*5d/)^^)    =   2((5d.5)f  -^)    =   2([5]N  -^) 


'AciTIBnbsuT    s,  zc^JBMi-tos   uioaj    'mon 


(a-9) 


I 


•     0    <    [~]Nt7    -    [5]N   ^ 


^p 


^^q^ 


saxiduiT    (lI-9)-b3   puB    'o    >   3   ^^^^   saxiduix   uaq:;    (6*9)*'^a 


(SI"9) 


^^    =   0   =    ^|M    c    ^   ^   0    =    ^ 

5e 


B^Bp   iBf^TUf    sq^    asooqo    puB 


'    0    >    [LilNgD   +    [IJ]M 


I 


^Bq:i 


This  is  an  ordinary  differential  inequality  for  N[|], 
subject  to  the  initial  conditions 

N[^]  =  N[n] 


|^N[§]    =   2aN[Ti].  (6.16) 


With   the   new  variable 


Eq.(6.15)    reduces   to 


,2 

2-^  >   0    ,  (6.18) 

dt 


and  the  initial  conditions  are 


y  =  0,  If  =  2a  .  (6.19) 


Integration  then  yields 


If  >  2a,   y  >  2at,   (for  t  >  0)  (6.20) 


or  equivalently  ^N[^]  >  2a  N[^], 


dt 

NEC]  >  N[n]  exp  2at   .  (6.21) 
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Using  these  tv;o  inequalities  in  Eq.  (6.1^1),  v;hlch,  in 
terms  of  the  kinetic  energy  K  =  N[3^/9t],  can  be  rewritten 
as 


K  >  (^  N[C])^/4N[C],  (6.22) 


•dt 


we  finally  obtain 


K(t)  >  K(0)  exp(2at)  (6.23) 


Hence 


lim  K  exp(-2at)  =  K(0)>0 

t-»-oo 


Sufficiency 

Taking  the  scaler  product  of  Eq.(6.6)  with  8u/8t  we  obtain 

N[8u/9t]  +  W[u]  =  C  (6.24) 

This  is  analogous  to  energy  conservation,  even  though  C  is 
not  the  energy  (it  does  not  even  have  the  dimension  of  an 
energy).   If  we  assume  that  W^  >_  0,  then 

W[u]+a2N[u]  >  0.  (6.25) 


HO 


After  eliminating  W  with  the  help  of  Eq.(6.2A)  this  inequality 
becomes 


2 
N[3u/3t]  <  C+o    N[u]   •  (6.26) 


On  the  other  hand,  the  Inequality  (6.14),  when  written 
for  u  rather  than  for  E,>    is 


(f^N[u])2  <  i4N[u]N[3u/3t]  .  (6.27) 


Combining  the  two  inequalities  (6.26)  and  (6.27),  and 
observing  that  N[u]  is  the  kinetic  energy,  we  obtain 

II  <      2>foV~T~CK  (6.28) 


Since  a^K^+CK  <_   ((jK+^)^,  this  implies 


II  <  2aK  +  C/a, 


and,  upon  integration. 


K  <  (K(0)+'^t)  exp(2at) 
~       a 


^1 


Kence,  for  a>a. 


11m  K  exp(-2at)  =  0 


For  cf  =  0 ,  Eq.(6.28)  implies  that  K  Is  bounded  by  a 
quadratic  function  of  time,  consistent  with  our  conjecture 
that  nonexponentlal  growth  of  K  Is  always  possible. 


il2 


7  .   Sufficient  Stability  Criteria 

The  advantage  of  the  energy  principle  lies  in  the  pos- 
sibility of  using  inequalities.   Therefore,  rather  than  solving 
the  Euler  equation  for  the  minimum  of  W^^  (which  would  be  as 
difficult  as  solving  the  equations  of  motion),  one  can  derive 
either  necessary  or  sufficient  conditions  for  stability. 
A  necessary  condition  is  optimistic  in  that  an  equilibrium 
may  be  unstable  even  if  the  condition  is  satisfied;  a  suffi- 
cient condition  is  pessimistic  in  that  an  equilibrium  may  be 
stable  even  if  the  condition  is  violated.   Correspondingly, 
necessary  conditions  can  be  used  only  to  shov/  that  an 
equilibrium  is  unstable,  and  sufficient  conditions  can  be 
used  only  to  show  that  an  equilibrium  is  stable. 

In  the  expression  (5-15)  for  W,  all  terms  are  positive 
except  the  one  involving  the  quantity  A  [Eq.(5.l6)].  Hence 
A  >_  0 ,  or  explicitly, 

(JxVp)- (B-V)Vp  <  0  (7.1) 

is  sufficient  for  exponential  stability.   In  order  to  exhibit 
stabilizing  equilibrium  properties,  v;e  apply  the  condition 
(7.1)  to  various  classes  of  equilibria. 

Let  us  first  assume  that  J.B  =  0.   Then,  from  Eq.(2.1) 

B^JxVp  =  -IvpI^b  , 
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and  the  condition  (7.1)  Is  equivalent  to 

B- (B-V)Vp  >  0  (7.2) 

The  Identity 

B- (B-V)Vp+Vp. (B.V)B  =  0  , 

which  follows  from  Eqs .  (2 .  S-'^ )  ,  nov;  Implies  that  the  condition 
(7.2)  Is  equivalent  to 

K  -Vp  <  0  ,  (7.3) 

where  1  is  the  curvature  of  a  magnetic  field  line, 

K  =  (b-V)b  ,   b  =  B/B  . 

Therefore,  an  equilibrium  is  exponentially  stable  If  no  current 
flov;s  along  the  magnetic  field  lines,  and  if  the  field  lines 
are  nowhere  concave  towards  Increasing  pressure.   Obviously, 
this  condition  is  never  satisfied  In  a  confined  plasma. 
Nevertheless,  It  gives  us  a  distinction  between  "good  curvature" 
and  "bad  curvature".   It  can  be  shown  that  an  equilibrium  Is 
stable  if  the  curvature  is  good  at  a  sufficiently  large  portion 
of  each  field  line.   The  stability  of  axlally  symmetric 
"multlpoles"  Is  a  consequence  of.  this.  ■  In  these  equilibria- 
the  pressure  has  its  maximum  at  the  separatrix  of  the  magnetic 
field.   Inside  the  separatrix,  the  curvature  is  good  everywhere; 


l\i\ 


I 


outside  the  separatrlx,  it  Is  good  near  the  x  point,  but 
bad  elsewhere.   It  can  be  shown  that  these  equilibria  are 
stable  as  long  as  the  plasma  does  not  extend  too  far  away 
from  the  separatrlx. 


SYMMETRY  AXIS 


SEPARATRIX 


WALL 


CONDUCTORS 
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The  curvature  criterion  (7-3)  is  still  satisfied  if 
the  curvature  is  identically  zero.   This  is  the  case  if  the 
magnetic  field  is  unidirectional, 

B  =  B(x,y)e^  . 

Then 

J  =  curl  B  =  VBxe 

~z 

Hence  J'B  =  0,  and  the  curvature  criterion  is  valid.   The 
pressure  balance  is 

V(p+|b2)  =  0     . 

An  equilibrium  of  this  kind  can  be  constructed  by  choosing 
some  domain  in  the  (x,y)  plane,  then  choosing  some  function 

B(x,y)  which  takes  a  constant  value  at  the  boundary,  and 

1  2 
finally  putting  p  =  -  pB  +const.   The  resulting  equilibrium 

is  exponentially  stable.   V/e  call  it  "confined"  if  p  has 

2 

its  minimum  at  the  boundary.   This  is  the  case  if  B  has  its 

maximum  at  the  boundary,  a  property  which  is  often  called 
"magnetic  well". 

Writing  the  pressure  balance  (2.1)  in  the  form 

Vp  =  (B-V)B-|VB^ 


or  equivalently 


B^K  +  (B-VB)b  =  V(p+1b^) , 


US 


we  have 


B^K-Vp  =  Vp-V(p+|b^)  , 


and  the  criterion  (7-3)  becomes 


Vp-VCp+ie^)  <  0  (7.^) 


This  implies  that 

Vp-VB^  <  0  (7.5) 

(unless  Vp  =  0),  thus  requiring  that  the  magnetic  field 

strength  decreases  In  the  direction  in  v;hich  the  pressure 

increases.   For  a  confined  plasma  in  which  the  pressure 

decreases  monotonically  from  its  maximum  to  the  wall,  the 

condition  (7-5)  Implies  a  magnetic  well. 

Good  curvature  and  magnetic  well  are  stabilizing 

properties  as  long  as  J-B  =  0  .   In  order  to  show  that 

this  need  not  be  so  in  general,  we  now  give  an  example  of 

a  whole  class  of  confined  equilibria  which  are  exponentially 

stable  independently  of  the  curvature  of  field  lines  and 

independently  of  the  gradient  of  the  magnetic  field.   These 

equilibria  are  characterized  by  8/9z  =  0  (plane  symmetry) 

and  J   =  0  (no  current  in  the  ignorable  direction).   Any 
z 

solenoldal  field  with  plane  symmetry  can  be  represented  as 


B  =  e^><ViJ;(x,y)+B  (x,y)e 


n 


Taking   the    curl   of   this,    one   finds 


J   =   e   Ai|j+VB   xe    , 


and  the  pressure  balance  becomes 


Ai|jViJj-B   VB   +[e    •(Vi|jxVB)]e      =    Vp 
z      z        z  z 


This  implies  that  both  B  and  p  are  constant  at  the  contours 
of  constant  4^,  and  that 


^*  =  i^(p4  Bz' 


Since  J   =  0  Is  equivalent  to 
z         ^ 

AiJ;  =  0,  (7.6) 

the  pressure  balance  in  these  equilibria  simply  reads 

p+^B^  =  const.  (7.7) 

2  z 

An  equilibrium  can  be  constructed  by  first  choosing  an 

arbitrary  domain  in  the  (Xjy)  plane,  then  solving  Eq.(7-6) 

with  the  boundary  condition  ti^  =  const,  and  finally  assigning 

values  of  both  p   and  B  ,  consistently  with  Eq.(7.7),  to 

z 

each  contour  of  constant  4).   In  order  that  the  solution  of 
Eq.(7.6)  is  not  just  the  trivial  one  (in  v/hlch  case  we  would 
deal  v;ith  a  unidirectional  magnetic  field  ^ 
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as  discussed  before)  we  have  to  choose  a  multiply  connected 
domain  with  different  values  of  ^   at  different  disjoint  parts 
of  the  boundary.   This  requires  current  carrying  Internal 
conductors  in  the  z  direction,  so  that  we  deal  with  a  straight 
multlpole.   The  stability  of  these  equilibria  follows  from  the 
fact  that  each  of  the  three  vectors  J,  Vp,  and  (B'V)Vp  is  in 
the  (x,y)  plane,  so  that  ( JxVp ) • (B- V) Vp  =  0. 

In  a  straight  circular  cylinder,  the  magnetic  field  is 
characterized  by  B   =  0,Bn(r),  B  (r),  where  (r,e,z)  are 
cylindrical  coordinates.   The  current  is  then  given  by 

J  =0,   J„  =  -—  B  ,  J   =  -  hS  r  B.   , 
r      6    drzz    rdr    6' 

and  the  pressure  balance  is 

The  sufficient  stability  condition  (7.1)  now  reduces  to 


J  B^  <  0 
z  6  — 


or 


^(rhl)<_0  (7.8) 


The  condition  (7-8)  cannot  be  satisfied  in  a  simply  connected 

2  2 
domain  unless  Bq  =  0 (unidirectional  field)  because  r  Bq  =  0 

2  2 

for  r=0,  but  r  B„  >  0  for  some  r  >  0. 
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Next  we  indicate  how  one  can  derive  better  (i.e.,  less  re- 
restrlctlve)  sufficient  stability  conditions.   From  Eq. 
(5.15), 


W  >  I  jJJ  dT(lF|^+Ax2) ,  (7.9) 


where 


JxVp 
F  =  Q+  ^= ^  X,  X  =  C-Vp 


Now,  from  Schwartz's  inequality. 


|Fp|Vpl  >  (F-Vp)' 


Further, 


F-Vp  =  Q-Vp  =  Vp'curl(C>=B)  =  div[  (^xB)  xVp]  =  div(XB)=  B-VX, 
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Hence 


|F|2  > 


(B-VX)' 

|vp|2 


and 


'M 


W  >  i  I  dT( 


(B-VX)' 


+  AX^). 


(7.10) 


Vp 


This  expression  has  two  advantages:   Firstly,  it  depends  only 
on  one  component  of  the  vector  trial  function  ^(viz.j  on  X); 
secondly,  only  derivatives  along  magnetic  field  lines  appear 
in  the  integrand. 

Consider  now  equilibria  either  with  an  ignorable  coordinate 
(like  with  axial  symmetry)  or  with  all  magnetic  field  lines 
closed  (constant  rational  rotation  number).   Then  all  physical 
quantities  are  periodic  along  magnetic  field  lines.   In  other 
words,  if  we  introduce  a  coordinate  s  increasing  along  the 
field  lines  according  to 


,     1  d£ 
^^  =  q-  -^  '   ^  = 


di 
B 


all  physical  quantities  are  periodic  in  s  with  period  unity, 
Equation  (7-10)  can  then  be  written  as 


W 


>-m< 


i-^ 


q  IVpT   ds 


^  I  +  AX^) 


X 


51 


where 


<•••> 


ds  ... 


Is  an  average  along  field  lines.   Obviously 


\a  Vp   ds      / 


.q  |Vp|   ds 

at  every  field  line  is  sufficient  for  W  >_  0 ,  and  hence  for 
exponential  stability.   The  condition  (7.11)  is  equivalent 
to  the  non-negativity  of  the  spectrum  of  an  eigenvalue 
problem  (as  the  non-negativity  of  the  MHD  spectrum  is 
equivalent  to  W  >  0).   This  eigenvalue  problem  is  given  by 


^    ^  X+AX  =  XgX  (7.12) 


ds   2|„  ,2  ds 
q  IVpl 

X(s+1)  =  X(s).  (7.13) 

Eq.  (7.12)  is  the  Euler  equation  of  the  functional  (7.11)  with 
the  normalisation 

dsgX^  =  1 

The  problem  (7.12-13)  is  considerably  easier  to  handle 
than  the  original  problem  not  only  in  that  it  involves  a 
second  order  ordinary  differential  equation  rather  than  a 
sixth  order  system  of  partial  differential  equations,  but 
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also  in  that  it  is  of  the  Sturmian  type:   the  spectrum  Is 
purely  discrete,  and  there  exists  a  lov;est  eigenvalue  X 
with  a  nonzero  eigenfunction.   A   >  0  at  every  field  line 
is  sufficient  for  stability.   To  determine  A   requires 
numerical  computations  in  general.   Unlike  for  the  original 
spectral  problem,  there  are  standard  procedures  for  this. 
It  is  also  possible  to  derive  lower  bounds  for 

A    analytically  (if  a  lower  bound  is  positive,  A   is  also 
positive,  and  the  equilibrium  is  stable).   The  condition 
A   >_  0  is  certainly  less  restrictive  than  the  previously 
discussed  condition  A>_0 .   In  order  to  get  a  rough  idea  of 
how  much  better  it  can  be  at  most  we  use  the  test  function 
X  =  1  in  Eq.(7.11)  to  obtain  <(a)  >_   0.   This  condition  is 
necessary  for  A  >_0 ,  and  hence  does  not  imply  stability. 
It  can  be  shown  that  \A/  <0  at  the  magnetic  axis  whenever 
J=fO  there.   Since  this  is  always  the  case  in  Tokamak-like 
equilibria,  the  present  sufficient  condition  is  not  good 
enough  for  investigating  the  stability  of  these  equilibria. 

Let  us  finally  derive  a  sufficient  condition  for  o 
stability.   From  Schwartz'  inequality, 

Ul^lvp!^  >  (c-vp)2  =  x^ 


Hence 


N  >  i 
--  ^ 


di 


X 


Vp 
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and,  using  Eq.(7.10),  we  find 


W 


2 


J  J  J     |Vp|  |Vpl 


For  equilibria  with  periodicity  along  field  lines, 

/lB:VXil  ^  ^X2^^2p  _xi.^N  >  0  (7.1^) 


at  every  field  line  is  sufficient  for  a  stability.   This  is 

2 

equivalent  to  a   >  -X  ,  where  A   is  the  lowest  eigenvalue 
^  —   o         o  ^ 

of  the  problem  (7.12-13)  with  the  weight  function 
g  =  p/|Vp|   .   In  other  words,  the  maximum  over  all  field 
lines  of  -X      is  an  upper  bound  for  the  maximum  growth  rate 
for  this  choice  of  the  weight  function. 

A  simpler,  but  weaker,  upper  bound  is  obtained  by 
dropping  the  first  term  in  Eq.(7-1^)-   This  yields 

2     „ 
(A+   ^  P^)X^   >  0 
|Vp|' 

as  a  sufficient  condition  for  a  stability.   This  inequality 
holds  for  all  X  if,  and  only  if. 


A+  ^-^  >  0  , 

.  |vp|'  - 


or 

^2  ,  _  A|7p|^ 
-       P 


Hence  the  maximum  of  the  r.h.s.  over  the  v/hole  domain  is  an 
upper  bound  for  the  maximum  growth  rate. 
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Sufficient  criteria  for  exponential  stability  have  never 
acquired  much  importance  because  they  appear  to  be  violated 
in  most  equilibria  of  practical  interest.   In  accordance  with 
this,  the  equilibria  which  we  have  shown  to  satisfy  a 
sufficient  criterion  either  are  not  confined  in  all  three 
spatial  directions  because  they  are  straight  or  have  internal 
conductors  (which  have  to  be  mechanically  supported  in  practice) 
Fortunately,  an  equilibrium  need  not  satisfy  a  sufficient 
criterion  in  order  to  be  acceptable  for  plasma  confinement, 
firstly  because  sufficient  criteria  are  overly  pessimistic, 
and  secondly  because  instabilities  with  small  grov^th  rates 
or  with  large  spatial  variations  are  not  always  dangerous. 
Therefore,  a  normal  mode  analysis  (i.e.,  an  analysis  of  the 
spectrum  together  with  the  eigenf unctions )  is  required 
to  indicate  whether  an  equilibrium  is  acceptable.   An 
intermediate  solution  of  the  problem  is  offered  by  evaluating 
sufficient  criteria  for  a  stability. 


?5 


8 .   Local  Stability 

Necessary  criteria  for  exponential  stability  can  be 
derived  by  considering  special  classes  of  trial  functions 
in  the  energy  principle.   If  a  necessary  criterion  is 
violated,  the  equilibrium  can  be  concluded  to  be  unstable; 
but  if  it  is  satisfied,  no  conclusion  can  be  drawn 
because  not  all  possible  trial  functions  were  considered 
in  its  derivation.   A  specific  necessary  criterion  arises 
if  the  subclass  consists  of  "localized"  trial  functions, 
i.e.,  functions  which  are  nonzero  only  in  an  arbitrarily 
small  neighborhood  of  a  pressure  surface.   The  equilibrium 
is  called  "locally  unstable"  if  this  criterion  is  violated. 
It  turns  out  that  this  criterion  takes  different  forms 
depending  on  the  topology  of  magnetic  field  lines.   More 
precisely,  the  version  which  is  valid  for  sheared 
equilibria  does  not  reduce  to  the  version  which  is  valid 
for  equilibria  with  all  field  lines  closed  in  the  limit  of 
vanishing  shear.   This  somewhat  paradoxical  fact  indicates 
that  the  threshold  for  exponential  stability  is  a  discon- 
tinuous surface  in  parameter  space,  and,  therefore,  that 
exponential  stability  need  not  always  be  a  very  meaningful 
concept  in  practice.   We  give  the  derivation  of  the  local 
stability  criteria  for  circular  cylinders.. 

These  equilibria  were  described  in  Section  7-   In 
order  to  simulate  the  toroidal  structure  of  confined 
equilibria  we  identify  any  two  planes  z  ==  const  with 
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distance  L,  thus  deallne  with  a  topological  torus.  For 
simplicity,  we  put  L  =  2-n .  The  rotation  nuiriber  is  then 
given  by 


y(r)  =  B  /rB  .  (8.1) 


We  first  consider  the  functional  W  .   Since  9  and 
z  are  ignorable  coordinates,  this  decomposes  according  to 


a  ~    ,  ^    o  ~kin-' 
k,m 


if  we  Fourier-decompose  the  trial  function  \i    according  to 


§  =  I   §km(^)  expi(me+kz)   , 
k,m 


?i  (^)  =  n,   (r)e  +v,  (r)en+v;i  (r)e 
-km   '   ^kJn   -r   km   ~6  ''^km   ~z 

v/ith  Integers  m  and  k. 

It  follovjs  that  W   can  be  made  negative  if,  and  only  if, 
it  can  be  made  negative  for  at  least  one  Fourier-component 
of  E, .      Therefore,  we  may  consider  one  Fourier-component  at  a 

time.   If  we  write  VJ^  explicitly  in  terms  of  the 
components  of  E,    ,    observing  that  derivatives  v/ith  respect  to 
0  and  z  are  replaced  by  im  and  Ik,  we  see  that  no  derivatives 
of  V  and  w  appear.   The  minimization  vjith  respect  to  v  and 
w  thus  can  be  performed  algebraically  (i.e.,  by  completing 
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squares).   Since  the  algebra,  though  straightforward,  is 
extremely  involved,  we  do  not  present  it  here,  but  merely 
write  down  the  result  (omitting  subscripts  k  and  m) : 


mm  W  ^ 

o 
v,w 


dr[a  l(ru)'|^+b„l  i^^] 


AC 
^0    rD   ' 


+r[2--2(— -  -  kBg)^]   , 
r 


2    2 

A  =  a  p+F 


C  =  a^p(B^+Yp)+YP^' 


^    M  2  ,  /m   .  ,  2  >  ^ 
D  =  a  p  +(—  +  k  )C 


(8.2) 


mB, 


k. 


p  =  — ^kB   =  mB  (Vi+^) 
r    z     z   m 


For  exponential  stability,  we  take  a=0.   Then 


A  =  F"  , 


C  =  YP^ 


D  =  (—  +  k  )yP^   , 
r 


58 


and 


a   =    ^  F 


o     2^,  2  2 
m  +k  r 


However,  the  limit  as  a->-0  of  b   depends  on  whether  F=0  or 


not.   If  F  ^  0, 


,2  _  ,^6^'   ^^^^e 


b   =  F  +2Bo(-^)  -   o   o  o 


kB      mB  , 

+r[2-^^-o(~^  -  kB^)] 


2^,  2  2'  r     ^6' 
m  +k  r 


(which  is  independent  of  y);  if  F  :;  0,  there  appears  the 
additional  term 


^k^BgYp 


(m^+k^r^)(B^+yp)   . 


In  the  latter  case,  a   =  0.   From  Eq.(8.1)  it  is  clear 
that  F=0  at  a  particular  surface  if  the  perturbation  is 
constant  along  the  field  lines  there.   If  k  and  m  are  not 
both  zero  this  is  possible  only  if  these  field  lines  are 
closed.   If  the  equilibrium  has  shear,  then  for  a  given 
pair  (k,m)  there  may  or  may  not  be  some  isolated  values  of 
r  at  which  F=0  ("mode  rational  surfaces").   If  the  shear 
is  zero,  then  P  either  is  identically  zero  or  has  no 
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zeroes  at  all,  depending  on  whether  k  and  m  are  such  that  the 

perturbation  Is  constant  along  field  lines  or  not  (If  y  Is 

irrational,  F  Is  never  zero,  but  can  be  made  arbitrarily  small 

by  choosing  k  and  m  accordingly). 

Let  us  now  consider  localized  test  functions  u,  i.e., 

we  pick  same  surface  r=r  ,  put  r=r  +£x  v;ith  a  small  parameter 
^  o         o 

e,  and  then  assume  that  u  is  a  smooth  function  of  x  (rather 
than  of  r),  and  that  u=0  outside  the  Interval  -l<x<l.   Since 
u '  =  —  g^  ,  the  first  term  in  Eq.(8.2)  is  large  compared  to 
the  second  term  in  general.   Since  this  first  term  is  non- 


negative,  we  cannot  make  W  negative  in  this  way  unless  a 
(and  hence  also  F)  is  zero  at  r  .   We  thus  pick  a  mode  rational 
surface  if  the  equilibrium  has  shear.   If  the  shear  is  zero, 
the  present  procedure  does  not  yield  a  stability  condition 
unless  all  field  lines  are  closed.   In  the  latter  case  we  have 
to  choose  k  and  m  such  that  F^O. 

First  considering  equilibria  with  shear,  we  expand  W   to 
lowest  order  in  e  to  obtain 


W   «  i  a"(r  ) 
o    2   o   o 


^^1^  5^1'  ^  \^-o'> 


dx  I  u 


If  this  can  be  made  negative,  the  exact  W   can  also  be  made 

^      '  o 

negative  by  making  e  sufficiently  small.   Let  us  first  minimize' 

the  first  term  keeping  the  second  term  fixed.   By  Schwartz's  inequal-' 


^  I   du I  2 


dx|ur  >  i 


du*    1 2 
dx   =-^  ux 
dx 
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Denoting  the  r.h.s.  by  |l|  ,  v;e  further  have  jl]   >_  (Re  I)' 
But 


Re 


I  =  |(I+I«)  =  |_ 


,    .du*  ,du  it  V    1 


^  dl'^l 


1 
2 


dx|ul 


Hence 


,  I  du I  2    1 


dxjul' 


and 


W  "c  (i  a"(r  )+b  (r  )) 
O     0   o   o    o   o 


dx|ui 


Now  using  F(r  )  =  0,  we  find  that 


i  a"(r  )  +  b  (r  )  -  p '+?■  B^(H1)2 
8   o   o     o   o    ^0   z^y 


r=r 


o 


This  expression  depends  on  k  and  m  only  through  r  . 

If  it  is  negative  for  any  pair  (k,m),  or  equivalently ,  at  any 

rational  surface  r  ,  W   can  be  made  negative,  and  the  equilibrium 

o    o  o      J 

is  unstable.   Since  the  rational  surfaces  are  dense,  it  is 
negative  at  some  rational  surface  if  it  is  negative  anywhere. 
Hence 


p^^|b2(H_|2>o 


(8.3) 


is  necessary  for  exponential  stability,  provided  y  |  0  . 
This  is  Suydam's  condition. 
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Next  we  consider  equilibria  v/lth  all  field  lines  closed 
As  already  noted  we  then  have  to  take  F  =  0.   We  thus  have 


W  = 
o 


dr  ho  lu| 


where  b   Includes  the  extra  term  depending  on  y*   Using 
FeO  (and  the  equilibrium  equations)  we  find 


,     ,     t    2yo   ^e,  2 
b  -  (p  -^  -f—   — )y 


Obviously, 


^2^     r  — 
B  +YP 


(8.4) 


is  necessary  for  stability  provided  y=|=0 

In  the  limit  of  zero  shear  (y'-^O),  Suydam's  criterion  (8,3) 
becomes  p'^0,  and  thus  is  more  restrictive  than  the  criterion 
(8.4).   Nevertheless,  it  cannot  be  concluded  that  the  criterion 
(8,4)  is  weaker  than  Suydam's  criterion,  and  therefore  redundant, 
because  Suydam's  criterion  is  not  valid  whenever  the  criterion 
(8,4)  is  valid.   The  correct  conclusion  is  that  the  stability 
threshold  is  a  discontinuous  function  of  y'  at  y'  =  0,   For 
Instance,  Suydam's  condition  implies  that  a  confined  plasma 
(p'<0)  is  locally  unstable  if  the  shear  is  too  small,  whereas 
the  criterion  (8.4)  implies  that  It  can  be  locally  stable  if 
the  shear  is  zero. 
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The  local  criteria  (8.3)  and  (S.'l)  have  been  generalized 
to  arbitrary  equilibria  with  shear  and  with  closed  field  lines, 
respectively -Due  to  the  absence  of  symmetries,  the  derivation 
of  these  general  local  criteria  is  extremely  involved,  and 
cannot  be  given  here.   Since  the  explicit  form  of  these  cri- 
teria is  not  very  intuitive  either  (i.e.,  stabilizing  and 
destabilizing  effects  are  not  directly  visible),  we  do  not 
even  write  them  down.   Hov/ever,  it  should  be  remarked  that 
a  large  body  of  the  literature  on  MHD  stability  is  concerned 
v;ith  the  evaluation  of  these  local  criteria  for  specific 
equilibria.   For  this  reason,  it  is  of  extreme  importance  to 
have  some  information  about  the  nature  of  the  instabilities 
that  occur  if  the  local  criterion  is  violated.   For  circular 
cylinders,  this  information  is  derived  in  Sec.  9. 
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9 .   Normal  Modes 

We  derive  properties  of  the  unstable  spectrum  In 
circular  cylinders.   We  recall  that  normal  modes  are  obtained 
by  making  the  functional  W^  stationery.   For  circular 
cylinders,  this  functional  is  given  by  Eq.(8.2).   Making 
it  stationery  leads  to  the  Euler  equation 


^  a  ^  ru  =  b  u  (9.1) 

dr   a  dr       o 


With  the  boundary  conditions  u=0  at  r=0 (regularity  of  the 
vector  K)    and  at  r=R  (perfectly  conducting  rigid  wall), 

Eq.(9.1)  constitutes  an  eigenvalue  problem  v/ith  eigenvalue 

2 
parameter  a  . 

This  problem  is  non-standard  not  only  in  that  the 

2 

coefficients  a   and  b   depend  on  a   in  a  complicated  manner, 
a      o 

but  also  in  that  the  relevant  function  space  is  not  simply 
that  of  square  integrable  functions  u (  r)  ,  but  rather  that 
defined  by  the  square-integrability  of  the  vector  ^  (the 
underlying  norm  involves  derivatives  of  u,  and  thus  is  a 
Sobolev  norm).   Therefore,  the  well-developed  theory  of 
the  spectrum  of  linear  second  order  ordinary  differential 
operators  is  not  applicable.   For  instance,  there  does  not 
seem  to  be  a  simple  way  of  showing  directly  from  Eq.(9.1) 
that  the  spectrum  is  real  (a  fact  which  we  know  to  be  true 
from  more  general  considerations). 


en 


An  obvious  way  to  find  eigenvalues  Is  the  "shooting 
method":   One  solves  Eq.(9.1)  as  an  initial  value  problem 

with  initial  conditions  at  r=0.   The  value  of  u  at  r=R 

2         2  2 

then  is  a  function  of  o  ,  say  f(a  ).   Those  values  of  a 

for  which  f=0  are  the  eigenvalues  of  the  problem.   The 

eigenvalues  are  isolated  as  long  as  the  function  f  is 

well-behaved,  and  this  is  the  case  as  long  as  Eq.(9.1) 

is  non-singular  J  i.e.  a  /b   =f  0.   This  is  the  case  for 

2 

0  >0.   Hence  the  unstable  spectrum  is  discrete,  and  con- 

tlnua  may  occur  only  at  the  stable  side.   In  particular, 
unstable  eigenvalues  cannot  accumulate  anyv/here  except 
at  the  origin.   However,  it  should  be  noted  that  this 
conclusion  holds  only  for  fixed  mode  numbers  k  and  m.   If 
we  form  the  entire  spectrum  by  taking  the  union  of  all 
the   subspectra  pertaining  to  fixed  k  and  m,  we  may  obtain 
unstable  accumulation  points. 

We  now  derive  a  monotonicity  property  of  the  unstable 

2 

eigenvalues:  a   is  a  monotonically  decreasing  function  of 

the  number  of  zeroes  of  the  associated  eigenfunction.   This 
number  is  called  "radial  mode  number",  and  we  denote  it  by 
n.   Since  monotonicity  in  n  of  the  eigenvalues  is  a  well-known 
property  of  Sturmian  problems  [(viz . ,  eigenvalue  problems  of 
the  form  (fu')  '+gu  =  Au)],  it  is  sometimes  referred  to  as  the 
"Sturmian  property",  even  though  the  problem  (9-1)  is  not  of 
the  Sturmian  type.   In  order  to  demonstrate  the  Sturmian 
property,  we  denote  the  solution  of  the  above  initial  value 
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2  2 

problem  by  u(r,o  ),  and  the  n-th  zero  of  this  by  r  (a  ). 

2 
Since  the  eigenvalues  are  given  by  r  (a  )=R,  it  suffices 

2       2 
to  show  that  r   is  a  monotonic  function  of  a      for  a  >0. 
n 

2 

If  we  assume  analyticity  of  the  function  r  (a  ),  this 

follows  from  the  reality  of  the  spectrum  through  the  follow- 
ing well-known  theorem:  If  an  analytical  function  f(Z) 
is  real  on  some  curve  in  the  complex  Z  plane,  and  if  f  has 
a  local  extremum  along  this  curve  at  some  point  Z  ,  then 
f  is  also  real  on  at  least  one  second  curve  through  the 
same  point  Z  .   This  theorem  is  easily  proved  by  considering 

the  power  series  of  f  about  the  point  Z  .   In  our  present 

2 

application,  the  function  r'  (a  )  takes  the  place  of  f(Z), 

and  the  curve  on  which  f  is  real  is  the  real  axis.   If 
we  assume  that  r   is  non-monotonic ,  it  follows  that  r 

must  have  a  local  extremum.   The  theorem  then  implies  that 

2 

r   is  real  for  complex  values  of  a    .      Hence  the  equation 

2  2 

r  (a  )=R  has  complex  solutions  a   for  real  R.   Since  this 

contradicts  the  fact  that  the  spectrum  is  real,  the  function 
r  (a  ;is  monotonic. 

The  Sturmian  property  has  several  practical  consequences 
Firstly,  it  implies  that  the  n=0  mode  v/hich  is  the  mode  best 
described  by  ideal  MHD  because  it  has  the  smallest  possible 
spatial  variation,  is  also  the  most  important  one  because  it 
has  the  largest  growth  rate.   Secondly,  it  implies  that  the 
growth  rates  decrease  if  the  wall  radius  is  decreased^  an 
effect  which  is  called  "wall  stabilization".   Thirdly,  it 
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implies  that  an  equilibrium  Is  o  stable  (for  the  given 

2 
mode  numbers  m  and  k)  if,  and  only  if,  r,  (a  )>_Rj  or 

equivalently ,  if  no  solution  of  the  Euler  equation  (9-1) 

has  more  than  one  zero  in  the  open  interval  0<r<R. 

We  now  ask  the  question  what  happens  as  we  let 
n^=°.   The  Sturmian  property  and  the  discreteness  of  the 
unstable  spectrum  imply  that  there  are  two  possibilities: 
Either  the  eigenvalues  accumulate  at  the  origin,  or  there 
are  only  finitely  many  unstable  eigenvalues  (in  which  case 
the  counting  breaks  down  as  the  origin  is  crossed  because 
Eq.(9.1)  then  is  singular). 

As  one  would  anticipate,  accumulation  takes  place  in 
equilibria  with  shear  if  Suydam ' s  criterion  is  violatea 

at  a  mode  rational  surface.   To  demonstrate  this  it  suffices 

2 

to  show  that  the  number  of  zeroes  of  the  solutions  u(r,a  ) 

2 

of  Eq.(9-1)  tends  to  infinity  as  o   tends  to  zero  (if 

Suydam's  criterion  is  violated).   Eq.(9-1)  is  singular 

2 
at  the  mode  rational  surface  r=r   if  a  =C.   Near  this 

o 

surface,  it  takes  the  form 

(x^u  ')  *  =  cu,  (9.2) 


where 


X  =  r-r 
o 


2b 

o 

^2   ,^    2 
d  a  /dr 
o 


) 


r=r 
o 
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The  behavior  of  the  solutions  of  Eq.(9.1)  near  the 
singularity  Is  given  by  the  solutions  of  Eq.(9.2).   The 
latter  are 


u  oc  exp(B£ogix[  )  ^^2) 

/jirr 


where 


3 


:±/c+i 


The  function  (9-3)  has  Infinitely  many  zeroes  near  x=0 
if,  and  only  If, 

c  +  ^  <  0 

This  is  the  case  if  Suydam's  condition  is  violated. 

In  closed  line  equilibria,  Eq.(9.1)  is  singular  for 

2 

o  =0  only  if  k  and  m  are  such  that  k+my^O  (resonant  modes). 

Hence  no  accumulation  of  unstable  eigenvalues  takes  place 
for  non-resonant  modes.   For  resonant  modes,  as  one  v;ould 
again  anticipate,  the  origin  is  an  accumulation  point  of 
unstable  eigenvalues  if,  and  only  if,  the  local  stability 

criterion  (8.4)  is  violated.   It  again  suffices  to  show  that 

2  2 

the  number  of  zeroes  of  u(r,a  )  becomes  large  as  a   becomes 

2 
small.   For  k+m.M  =  0,  and  small  a  ,  Eq.(9.1)  becomes 


m  +k  r 
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since  the  coefficient  of  u"  is  small,  (and  positive)  the 

solutions  oscillate  rapidly  (and  thus  have  many  zeroes) 

wherever  b  <0.   This  Is  the  case  wherever  the  local 
o 

stability  condition  (8.^)  is  violated. 

How  can  we  now  understand  the  fact  that  there  are  two  distinct 
local  stability  criteria  for  u'=t=0  and  for  ii'=0?   This  fact 
seems  even  more  paradoxical  if  we  observe  that,  by  the  same 
argument  by  which  we  concluded  that  the  unstable  spectrum 
is  discrete,  each  unstable  eigenvalue  depends  continuously 
on  equilibrium  parameters.   The  answer  is  that  this  continuity 
is  non-uniform  in  the  mode  number  n.   Thus,  the  limits 

n->oo  and  m'-»-0  are  not  interchangeable:   Taking  first  n-»-o°, 

2 

and  then  u'-*0,  we  find  w  ^— 0  if  p'<0,  but  taking  first 

2 
y'-^O,  and  then  n-*-=° ,    vie    find  w  ->-0  if  the  criterion  (8.^) 

is  violated.   Similar  non-uniformities  can  be  found  for 

other  limits  (e.g.,  m->=o ,  k-^^ ,  |  p '  | -^o° ,  etc.). 

The  following  statements  suimnarize  what  we  have  learned 

about  the  unstable  spectrum: 

1)  The  unstable  spectrum  consists  of  point  eigenvalues, 
which  depend  continuously  on  the  equilibrium. 

2)  Each  unstable  point  eigenvalues  is  uniquely  characterized 
by  the  three  mode  numbers  k,m,  and  n. 

3)  For  fixed  k  and  m,  the  growth  rate  decreases  as  n 
Increases,  and  tends  to  zero  as  n->=°(if  the  n->«>  modes 
are  unstable  at  all). 
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^)      For  fixed  k  and  m ,  the  stability  of  the  n->-°°  modes 
depends  on  the  local  stability  criterion.   If  y'=fO, 
these  modes  are  unstable  if,  and  only  if,  Suydam's 
criterion  is  violated  at  a  mode-rational  surface; 
if  y  is  a  rational  constant,  they  are  unstable  if, 
and  only  if,  the  criterion  (S.'l)  is  violated  any- 
where in  the  plasma. 
From  a  practical  point  of  viev;,  only  those  unstable  modes 
are  important  whose  mode  numbers  are  small  ("gross  modes"). 
The  m=±l,  k=±l,  n=0  modes,  for  Instance,  are  of  particular 
importance  in  equilibria  with  shear  because  they  are 
frequently  observed  in  experiments.   From  the  above  statements 
it  is  clear  that  these  modes  are  unstable  if  Suydam's  criterion 
is  violated  at  a  y=±l  surface.   More  generally,  the  n=0 
mode  is  unstable  for  any  given  k  and  m  if  Suydam's  criterion 
is  violated  at  a  mode  rational  surface.   Suydam's  criterion 
thus  is  a  condition  for  gross  stability,  even  though  it  is 
derived  by  considering  localized  test  functions,  but  only 
at  the  low-order  rational  surfaces.   If  there  are  no  such 
surfaces,  as  it  is  the  case  in  equilibria  with  small  shear, 
Suydam's  criterion  is  not  related  to  gross  modes.   As  to  the 
relevance  to  gross  modes  of  the  criterion  (8.4)  for  closed 
line- equilibria,  there  is- no  restriction  to  particular 
surfaces  because  all  surfaces  are  mode  rational  if  k/m+y=0. 
Hence  the  criterion  is  relevant  as  long  as  y  is  a  low  order 
rational  number.   By  continuity,  it  is  still  relevant  if  y 
is  close  to  a  low  order  rational  number  (e.g.,  small  y  and 
low  shear)   Needless  to  say  that  these  practical  conclusions 
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should  also  apply  to  the  generalizations  to  arbitrary 
equilibria  of  the  local  criteria  (8.3)  and  (8.^). 

As  a  final  remark  we  stress  that  gross  MHD  stability 
of  experimental  plasmas  is  a  quantitative  numerical  problem. 
The  local  criteria  are  only  necessary,  but  not  sufficient 
for  gross  stability,  and  thus  yield  only  a  bound  for  the 
threshold.   The  actual  threshold,  unlike  the  local  criteria, 
is  not  related  to  a  singularity  of  the  equations,  and  thus 
cannot  be  determined  analytically.   This  is  true  even  more 
for  the  growth  rates  of  gross  modes.   And  if  there  is  any- 
thing in  the  linearized  MHD  theory  which  is  relevant  to 
plasma  confinement,  then  these  are  the  growth  rates  of  gross 
modes . 
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This  report  was  prepared  as  an  account  of 
Government  sponsored  work.   Neither  the 
United  States,  nor  the  Administration, 
nor  any  person  acting  on  behalf  of  the 
Administration : 

A.  Makes  any  warranty  or  representation, 
express  or  implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  information  contained  in  this  report, 
or  that  the  use  of  any  information, 
apparatus,  method,  or  process  disclosed 
in  this  report  may  not  infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
the  use  of,  or  for  damages  resulting  from 
the  use  of  any  information,  apparatus, 
method,  or  process  disclosed  in  this 
report . 

As  used  in  the  above,  "person  acting  on  behalf 
of  the  Administration"  includes  any  employee 
or  contractor  of  the  Administration,  or 
employee  of  such  contractor,  to  the  extent 
that  such  employee  or  contractor  of  the 
Administration,  or  employee  of  such  contractor 
prepares,  disseminates,  or  provides  access  to, 
any  information  pursuant  to  his  employment  or 
contract  with  the  Administration,  or  his 
employment  with  such  contractor. 
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